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Abstract: We describe the impact of the full one-loop EW terms of 0{asOi\y^ entering 
the electron-positron into three-jet cross-section from ^/s = Mz to TeV scale energies. We 
include both factorisable and non-factorisable virtual corrections, photon bremsstrahlung 
but not the real emission of and Z bosons. Their importance for the measurement of 
as from jet rates and shape variables is explained. 
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1. Introduction 

Strong (QCD) and Electro- Weak (EW) interactions are two fundamental forces of Nature, 
the latter in turn unifying Electro-Magnetic (EM) and Weak (W) interactions in the Stan- 
dard Model (SM). A clear hierarchy exists between the strength of these two interactions at 
the energy scales probed by past and present high energy particle accelerators (e.g., LEP, 
SLC, HERA and Tevatron): QCD forces are stronger than EW ones. This is quantitatively 
manifest if one recalls that the value of the QCD coupling, as, measured at these machines 
is much larger than the EM one, oeMj typically, by an order of magnitude. This argument, 
however, is only valid in lowest order in perturbation theory. 

A peculiar feature distinguishing QCD and EW effects in higher orders is that the 

latter are enhanced by (Sudakov) double logarithmic factors, ln^(-7^), which, unlike in 

w 

the former, do not cancel for 'infrared-safe' observables ^, ^. The origin of these 

'double logs' is well understood. It is due to a lack of the Kinoshita-Lee-Nauenberg (KLN) 
||5| type cancellations of Infra-Red (IR) - both soft and collinear - virtual and real emission 
in higher order contributions originating from (and, possibly, Z: see Footnote 1) 
exchange. This is in turn a consequence of the violation of the Bloch-Nordsieck theorem 
Q in non-Abelian theories The problem is in principle present also in QCD. In practice, 
however, it has no observable consequences, because of the final averaging of the colour 
degrees of freedom of partons, forced by their confinement into colourless hadrons. This 
does not occur in the EW case, where the initial state has a non-Abelian charge, dictated 
by the given collider beam configuration, such as in e~^e~ collisions. 

These logarithmic corrections are finite (unlike in QCD), as the masses of the weak 
gauge bosons provide a physical cut-off for and Z emission. Hence, for typical ex- 
perimental resolutions, softly and collinearly emitted weak bosons need not be included in 
the production cross-section and one can restrict oneself to the calculation of weak effects 
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originating from virtual corrections and affecting a purely hadronic final state . Besides, 
these contributions can be isolated in a gauge-invariant manner from EM effects at 
least in some specific cases, and therefore may or may not be included in the calculation, 
depending on the observable being studied. As for purely EM effects, since the (infinite) 
IR real photon emission cannot be resolved experimentally, this ought to be combined with 
the (also infinite) virtual one, through the same order, to recover a finite result, which is 
however not doubly logarithmically enhanced (as QED is an Abelian theory). 

In view of all this, it becomes of crucial importance to assess the quantitative relevance 
of such EW corrections affecting, in particular, key QCD processes studied at past, present 
and future colliders, such as e~^e~ — > 3 jets. 



2. Three-jet Events at Leptonic Colliders 

It is the aim of our paper to report on the computation of the full one-loop EW effects 
entering three-jet production in electron-positron annihilation at any collider energy via the 



subprocesses e~^e~ —>■ Z —>■ qqg^. Ref. |10| tackled part of these, in fact, restricted to the 
case of and Z (but not 7) exchange and when the higher order effects arise only from 
initial or final state interactions (these represent the so-called 'factorisable' corrections, 
i.e., those involving loops not connecting the initial leptons to the final quarks or gluons, 
which are the dominant ones at ^/s = Mz, where the width of the Z resonance provides 
a natural cut-off for off-shellness effects). The remainder, 'non- factorisable' corrections, 
while being typically small at ^/s = Mz, are expected to play a quantitatively relevant 
role as ^/s grows larger. Since, here, we study the full set the one-loop EW corrections, 
we improve on the results of Ref. [I0| in two respects: (i) we include now all the non- 
factorisable terms; (ii) we also incorporate previously neglected genuine QED corrections, 
including photon bremsstrahlung. In contrast, we refrain here from computing and 
Z boson bremsstrahlung (just like in ||l^), as we will argue that this may not enter the 
experimental jet samples. (For a study of the impact of the emission of real massive gauge 



bosons in a variety of high energy processes, see Ref. |11|, albeit in hadronic collision.) 

Combining the aforementioned logarithmic enhancement associated with the genuinely 
weak component of the EW corrections to the fact that as steadily decreases with energy, 
unlike oeWi in general, one expects one- loop EW effects to become comparable to QCD 
ones at future Linear Colliders (LCs) running at TeV energy scales, like those available 
at an International Linear Collider (ILC) or the Compact Linear Collider (CLIC)^. In 
contrast, at the Z mass peak, where logarithmic enhancements are not effective, one-loop 
EW corrections are expected to appear at the percent level, hence being of limited relevance 



""^By doing so, double logarithms also arise from virtual Z exchange, despite there being no flavour 
change, simply because KLN cancellations are spoilt by the remotion of real Z emission. 

^See Ref. Q for the corresponding one-loop corrections to the Born process e^e~ —> qq and Ref. for 
the ~ ni component of those to e^e~ — > qqgg (where rif represents the number of light flavours). 

^For example, at one-loop level, in the case of the inclusive cross-section of e^e~ into hadrons, the QCD 
corrections are of Of — ), whereas the EW ones are of Of "f^ In'^ -r^), where s is the collider centre-of-mass 
(CM) energy squared, so that at ^/s « 1.5 TeV the former are identical to the latter, of order 9% or so. 
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(d) 



(f) 



Figure 1: Internal and external bosonic self-energy graphs. The shaded blob on the wavy lines 
represents all the contributions to the gauge boson self-energy and is dependent on the Higgs mass, 
as detailed in the graphical equation. These corrections apply to the photon self-energy, the Z 
boson self-energy and the photon-Z mixing. The gauge bosons inside the loops are and the 
scalars are charged Goldstone bosons. In the case of the self-energy of the Z, in graphs (d)-(f), the 
gauge boson inside the loop can be a Z and the scalar a Higgs boson. 



at LEPl and SLC, where the final error on as is of the same order or larger [|^, but of 
crucial importance at a GigaZ stage of a future LC |1C], where the relative accuracy of as 



measurements is expected to be at the 0.1% level or better |14]. On the subject of higher 
order QCD effects, it should be mentioned here that a great deal of effort has recently been 



devoted to evaluate two- loop contributions to the three-jet process |15| while the one-loop 



QCD results have been known for quite some time |16]. 

As intimated, in the case of e~^e~ annihilations, the most important QCD quantity to 
be extracted from multi-jet events is as- The confrontation of the measured value of the 
strong coupling constant with that predicted by the theory through the renormalisation 
group evolution is an important test of the SM. Alternatively, it may be an indication 
of new physics, when its typical mass scale is larger than the collider energy, so that the 
new particles cannot be produced as 'real' detectable states but may manifest themselves 
through 'virtual' effects. Not only jet rates, but also jet shape observables would be affected. 

The calculation we are presenting here involves the full one-loop EW corrections to 
three-jet observables in electron-positron annihilations generated via the interference of 
the graphs in Figs. with the tree-level ones for e^e~ — * 7*,^ — > qqg- Hence, our 
calculation not only accounts for the aforementioned double logarithms, but also all single 
ones as well as the finite terms arising through the complete ©(asag^y)- We will account 
for all possible flavours of (anti)quarks in the final state, with the exception of the top 
quark. The latter however appears in some of the loops whenever a bbg final state is 
considered. 

We will show that all such corrections can range from a few percent to a few tens of 
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Figure 2: Internal and external fermionic self-energy graphs. The gauge boson mass, M, can take 
the values A, Mz or Mw- For the external self-energies, we display only the self-energy insertion 
on the incoming positron (c), (d) or outgoing antiquark (e), (f). There is an identical contribution 
from self-energy insertions on the electron or quark, but each carries a combinatorial factor of i. 

percent, both at = Mz and LC energies, depending on the observable under study, with 
the QED component being preponderant at low energy and the weak one emerging more 
and more as the latter increases. Altogether, while their impact is not dramatic in the 
context of LEPl and SLC physics at a GigaZ stage of future LCs they ought to be taken 
into account in the experimental fits. This is even more the case of future LCs running at 
and beyond the TeV range^. 

The plan of the rest of the paper is as follows. In the next Section, we describe the 
calculation. Then, in Sect. ^ we present some numerical results (a preliminary account of 
which was given in Ref. [||]). We conclude in Sect. |. 

^Another relevant phenomenological aspect of our calculation would pertain to the case in which the 
bremsstrahlung photon is resolved, so that a three-jet plus one-photon sample can be defined and compared 
to a four-jet sample in order to extract genuine non-Abelian QCD effects, as detailed in - and references 
therein - but we postpone this analysis to a future work. 
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(a) (b) 




(e) (f) 



Figure 3: Vertex graphs. Corrections on the (anti)quarli line. In graphs (a) and (b) the gauge 
boson mass inside the loop must be Mw whereas for graphs (c) to (f) the gauge boson mass M can 
take the values A, Mz or Mw- 

3. Calculation 

In Ref. [^], which concentrated only on the factorisable corrections to three-jet production 
(i.e., neglecting any interaction between the incoming electron-positron pair and the final- 
state (anti)quarks), the method of helicity Matrix Elements (MEs) was adopted. In the 
present case such an approach is not convenient because we have to address the problem 
of IR divergent contributions from loops in which a photon is exchanged. As intimated 
already, this means that we also need to consider the real photon (bremsstrahlung) contri- 
butions, summed over the helicities of the emitted photon. Moreover, since in our numerical 
simulations, as we will see later, we regularise the collinear singularities with finite fermion 
masses, the introduction of these terms would complicate the expressions obtained within 
the helicity formalism. For the virtual corrections, we therefore adopt the more tradi- 
tional approach of considering all possible interferences between the one-loop and tree-level 
graphs^. We have organized the former into "prototype graphs", shown in Figs. which 
^We have however verified that, limitedly to the case of factorisable corrections, we can reproduce the 
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(c) (d) 

Figure 4: Vertex graphs. Corrections on tlie electron/positron line. In graphs (a) and (b) the 
gauge boson mass mass M can take the values A, Mz or Mw whereas in graphs (c) and (d) inside 
the loop there must be Mw- 



were calculated as general functions of couplings and internal masses, which appear inside 
Passarino-Veltman (PV) functions [^^, which are then determined numerically with the 
appropriate couplings and internal masses. In doing so, we have made extensive use of the 
program. 



FORM 21 



In anticipation of an electron-positron collider in which the incoming beams can be 
polarised, we have however inserted a helicity projection operator into the electron line 
and obtained separate results for left-handed and right-handed incoming electrons^. For 
genuinely weak interaction corrections, this is of particular interest, since such corrections 
violate parity conservation. Unfortunately, this occurs already at tree-level, owing to the 
contribution from exchange of a Z boson, but these higher order corrections are also pe- 
culiarly dependent on the incoming lepton helicity and thus one would expect the two 
parity-violating effects be distinguishable after the collection of sufficient events. (We will 
devote a separate paper to the study of polarised incoming beams.) 

In this connection, it is worth mentioning that as we are dealing with weak interactions, 
which involve axial couplings, and as we furthermore wish to distinguish between left- and 
right-helicity incoming electrons, the only unambiguous way to handle loop integrals is 
to use the 'modified' Dimensional Reduction (DR) as opposed to the 'modified' Minimal 
Subtraction (MS) renormalisation prescription. (We use ^ = Mz as the subtraction scale 
of ultraviolet divergences.) It is therefore necessary to determine the value of the fine- 



(unpolarised) results of Ref. . 

® As we are taking massless incoming fermions, the helicity of the positron is simply the opposite to that 
of the electron. 
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(f) (g) 

Figure 5: Box graphs. In the non-factorising corrections (d)-(f), the gauge boson masses Afi, M2 
can each take the values A or Mz or they can both be Mw- 

structure constant oem in that scheme at the subtraction point ^ = Mz^ ■ For the light 
fermion (including the 6-quark) contributions we can integrate the /3-function through 

^ We emphasise that we are not using the complex mass subtraction scheme of Ref. , but rather the 
DR one in which all counterterms are taken to be real. For the internal gauge boson self-energies the mass 
subtraction is the on-shell self-energy. Provided this is taken in conjunction with the corresponding values 
for the EM coupling, qem, and the weak mixing angle, Ow (as discussed here), this is an equally consistent 
scheme in which the imaginary parts of the renormalised couplings in the scheme of are reproduced by 
the imaginary parts of the subtracted gauge boson self-energies. 
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Figure 6: Pentagon graphs. The gauge boson masses can take the same values as in the case of 
the non-factorising box graphs. 



the thresholds at = 4m?, for all charged fermions with mass 2mi < Mz- There is 
a difficulty here with the light quarks as one expects substantial QCD corrections and 
the mass thresholds are well below the regime where perturbative QCD is reliable. This 
problem is usually addressed by considering the total cross-section for electron-positron 
annihilation in the resonant region and applying a finite energy sum rule. This has been 
carried out in detail in Ref. |23]. The upshot of this numerical analysis is that the effect of 



all light fermions can be simulated by using ordinary perturbation theory with a threshold 
of 15 MeV for the u- and d-quarks and 1 GeV for the s-quark, yielding a value aEM^(-^l) ~ 
1/128.2, in the MOMentum (MOM) subtraction renormalisation scheme. However, there is 
a finite difference between this prescription and a genuine DR determination arising from: 



1, 



contributions to the photon self-energy from loops of gauge bosons (and their 
attendant Goldstone bosons); 



2. a contribution from the non- vanishing Z photon mixing part of the self-energy at 
zero momentum; 

3. a contribution from the loop correction to the EM vertex involving internal VF^'s, 
which does not cancel by virtue of the QED Ward identity, owing to the non-Abelian 
nature of the coupling of the photon to charged gauge bosons. 



These corrections have been considered in Ref. |24| in the MS scheme. In the DR scheme 
we find that these give a total contribution of 

A("iS(M|)) =-^^iMM),„(^), ,3.) 

where sw and cw are sm9]y and cos 6w, respectively. This introduces a negligible correc- 
tion. 

Finally, we need to account for the difference between the MOM and the DR schemes, 
from the light fermions (including the 6-quark) , and a contribution from the t-quark, which 
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lead to the more substantial corrections 



EM 



(M; 



Svr 



(3.2) 



where the sum is over all fermions whose electric charge is gj e, except the i-quark. This 
gives us a final value of 



DR 

EM 



{Ml 



1 



(3.3) 



127.7 

Loops containing one or two photons give rise to IR divergences. For all of these 
divergent loop integrals we have a means of isolating the IR part in the DR scheme in 
4 — 2e dimensions. This is explained in detail in the appendix. For all of these IR divergent 
integrals we have checked numerically that the introduction of a common mass, A, for both 
the internal fermion lines and the internal photon lines reproduces these integrals upon the 
replacements 

7r-T(e) ^ ln(A2), 



TT 



ln^(A^ 



e 2 

provided A is chosen sufficiently small, but not so small as to introduce numerical insta- 
bilities. In each case we have identified a substantial plateau region in A where the two 
methods of calculation agree to a very high degree of accuracy. 

Having established that this works, we then abandon the use of dimensional regularisa- 
tion for the case of IR divergences and simply insert a mass regulator A for the photon and 
an independent mass mj for all fermions. This is also done in the case of the bremsstrahlung 
contribution before integrating over the phase space for the emitted photon. The reason 
why the mass regularisation of IR (i.e., soft/collinear) singularities works is linked to the 
fact that the tree-level diagrams only involve neutral currents, so that the EM and purely 
weak corrections are separately gauge invariant, as mass regularisation is well known to 
work in QED. By varying both of these parameters and noting the insensitivity of the final 
results to these changes, we have checked the expected cancellations of these divergences 
between the virtual correction and real emission. 

A new feature of this calculation, which was not present in the case of factorisable 
corrections only dealt with in ||l^, is the occurrence of pentagon graphs, as shown in 
Fig. |6|, which arise from interactions with two gauge bosons exchanged in the s-channel 
and the gluon emitted from an internal (anti)quark. Such graphs involve five-point PV 
functions with up to three powers of momenta in the numerator. We have handled these in 
two separate ways (with two independently developed codes), in order to check for possible 
numerical instabilities. In the first case the integrals are simply evaluated using routines in 
LoopTools v2.2 [E^], for which these tensor-type pentagon integrals are now implemented 



according to the reduction formalism of Ref. |26|*. In the other we use the standard PV 
reduction to express scalar products of the loop momentum k with external momenta pi 



*We implemented also directly in an independent FORTRAN routine the expressions for the five-point 
functions of Ref. |26|, finding agreement with LoopTools up to available digits precision. 
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in terms of the denominators of propagators adjacent to the i-th vertex. This reduction 
is carried out exhaustively until only scalar pentagon integrals appear. The latter are 



available in the library FFl .9 |27]. A comparison of the numerical results provided by the 
two codes for a sample of phase space points yielded satisfactory agreement between the 
two methods. 

The square MEs and the interferences for the bremsstrahlung process have been cal- 
culated directly using ALPHA pM] and checked against both MadGraph |2£] and the results 



based on the codes of Ref. |3C], finding perfect agreement between all these implementa- 
tions. 

Regarding the integration over the photon phase space, this has been split into two 
regions 

1. < AE. 

2. E^ > AE. 

For the first region, the MEs have been approximated using the eikonal approximation 
in which a fermion with momentum p emitting a photon (with polarisation e) has an 
associated vertex 2p ■ e and the photon momentum is neglected in all numerators. Here, 
the phase space integration has been performed analytically assuming a photon mass A. In 
the second region, the full MEs are employed but the photon mass is set to zero. Here, we 
integrate over all the phase space available to the photon. The integrations for both the 
three- and four-particle final states were performed numerically with Monte Carlo (MC) 
techniques, so that any desired cut on the final state can be imposed. 
The three- and four-particle Lorentz invariant phase space reads as 

^ = W^^^^ ^ (^^^ -p,-p,-p,l 

^ - 16(2^^^^^ ^ ^P^^ -P,-P,-P,-P,), (3.4) 

respectively. We choose as five (eight) independent integration variables for the three (four) 
body phase space the angles of the quark, the energy and the angles of the gluon (and the 
energy and the angles of the photon). With this choice, Eq. ( p.4| ) can be written as 

8(2^)5 |ii;,- + i?,(i + Pg-p,/|p,|2) I 3 ' ^' 
d8$4 = — i-— — ———dE.d^nadE^d^n^d^n,. (3.5) 

16(27r)8|ii;,- + i?,[l + (p, + p^).p,/|p,|2]| 9 9 y J , y J 



In Eq. (^]^), all the kinematical quantities are derived from the independent variables by 
imposing four-momentum conservation and mass-shell relations and are calculated in the 
rest frame of the incoming e^e~ system. 

In order to perform the numerical integration over the phase space, the independent 
variables are sampled according to the peaking structure of the integrand function. Stan- 
dard multi-channel importance sampling techniques are used. Concerning the importance 
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A2 (GeV2) 


5 • 10-10 


1 • 10-13 


(yv + f^s (pb), v^ = 90 GeV 
(TV + (Ts (pb), ^/s = 350 GeV 


-4181.6 ±5 
-1.2293 ±0.005 


-4181.2 ±5 
-1.2293 ±0.005 



Table 1: Variation of the soft plus virtual cross section ay + cfs for two different values of A for 
= 90 and 350 GeV. Here, only order oem QED contributions are included, the cuts are specified 
in the text and 2^Ej^/s is lO-''. The sum over the final-state quark fiavours is taken. 





10-3 


10-^ 


(Treai (pb), ^ = 90 GeV 


5054. ± 1. 


5050. ± 2 


o"reai (pb), ^fs = 350 GeV 


1.978 ±0.001 


1.979 ±0.001 



Table 2: Variation of the four-body cross section (iroai for two different values of Ai?, wherein A'^ 
is fixed to the value 10- ^'^ GeV^. Again, only order q;em QED contributions are included, ^fs is 90 
and 350 GeV, the cuts are specified in the text and the sum over the final-state quark flavours is 
taken. 



sampling of the three- and four-body phase space collinear configurations, our treatment 
closely follows the multi-channel approach described in appendix A. 3 of Ref. |31]. 

As intimated, we have checked that the final result is insensitive to variations of (small 
values) of not only the fictitious photon mass A but also the separator Ai?. Table |l| shows 
the sum of the virtual, cry, and soft, as-, cross section (integrated on from A to AE'), 
while spanning a wide range of values for A. Table |2| shows the real radiation cross section, 
Creai) where the integration on is split in the soft part A < E'^ < AE' and the hard 
part A£' < E^ < E^^^. (Parameters and cuts on the hadronic system are defined in the 
following Section.) 

When the real photon emission process is combined with the virtual corrections we 
have checked that the limit A — > may be taken, whilst keeping the fermion mass rrij 
small but non-zero. The collinear divergence obtained when — > cancels with the 
well-known exception of an overall large logarithm (ln(s/mj)), which is associated with 
the Initial State Radiation (ISR) induced by the incoming electrons and positrons. 

It is well known that, in the case of QCD, such a remnant collinear divergence is 
absorbed into the momentum dependence of the Parton Distribution Functions (PDFs). 
In the case of electron-positron colliders this large correction is always present, but it is 
universal to all processes and would therefore tend to mask the purely weak corrections. 
However, for sensible numerical results, it has to be accounted for to all orders of perturba- 
tion theory, e.g., within the so-called electron/positron structure function formalism |^2|| , 



which automatically resums in QED all Leading Logarithmic (LL) terms. In Ref. |33] a 
method of combining consistently resummed LL calculations with exact ©(oem) ones has 
been devised both in additive and factorisable form. Here, we will adopt the additive 
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approach, which amounts to write the differential cross section as 



da = do-LL - dal^ + dcJexact > (3-6) 
where dai,h and da^^^ are defined by 

dfiLL = J dxidx2D{xi, s)D{x2, s)daQ{xiX2s), 

da^i^ = J dxidx2[D{xi, s)D{x2, s)]adaQ{xiX2s) (3-7) 

and where D{x,Q'^) is the electron structure function whilst [Z)(xi, s)D(x2, s)]a is the 
order oem term of D{xi, s)D{x2, s). In Eq. ( |3.6D , the order oem part of the up-to-all-order 
resummed cross section in LL approximation is replaced by the exact order oem correction 
or, in other words, the higher-order (beyond order c^em) LL corrections are added to the 
exact order oem calculation. 

There remains one kinematic region where our results would not be numerically reli- 
able. In the case of real photon emission, the aforementioned ISR can lead to the situation 
in which the remaining sub-energy of the incoming electron-positron pair is close the the 
Z mass and hence close to the pole of the propagator in the case of Z boson s-channel 
exchange. In order to avoid this the width, Tz, of the Z boson has been included in the 
propagator. For consistency, this means that the same width has to be included in the 
Z propagator for the virtual corrections. This can be done in a consistent way in the 
DR scheme used in this paper, in which renormalised couplings remain real, as well as 



the the complex mass scheme [22], where the W and Z masses are defined as the loca- 
tions of the propagator poles in the complex plane, and the couplings become complex. 
The essential ingredient for the evaluation of virtual corrections is the ability to compute 
one-loop integrals with complex internal masses. The package LoopTools v2.2 imple- 
ments complex masses only for two-point, three-point and IR-singular four-point functions. 
We implemented in LoopTools the general expression for the scalar four-point function of 
Ref. Q], valid also for complex masses. Particular attention has been devoted to the 
occurrence of numerical instabilities in certain regions of phase space because of strong 
cancellations. The new routine for the scalar four-point function with complex masses 
has been tested with extended coinplex*32 precision, finding perfect numerical agreement. 
Actually, for the present study, the implementation of the complex mass scheme is only 
partial since we have kept real couplings and thus our scheme is equivalent to the fixed 



width scheme |35| implemented at the full one-loop level. In principle this could produce 
a bad high energy behaviour of the corrections. We tested the stability of our predictions 
at high energy by switching off the W/Z width in the virtual and real corrections, check- 
ing that the effect of the width is smaller and smaller as the energy raises up to 1 TeV: 
the differences on the pure virtual corrections with and without W/Z width (and setting 
A2 = 5 • 10-10 GeV2) are 0.5% at =350 GeV, 0.2% at ^/s =700 GeV and 0.1% at 
=1 TeV, in units of the Born cross section. 

We have neglected the masses of light quarks throughout. However, in the case in 
which the final state contains a bb pair, whenever there is a boson in the virtual loops. 
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account had to be taken of the mass of the top (anti)quark. Furthermore, in such cases, 
it was also necessary to supplement the graphs shown in Figs. |^-^ with the corresponding 
diagrams in which one or two internal bosons are replaced by charged Goldstone 
bosons, as appropriate. Having done all this, we are therefore in a position to present the 
results for such '6-jets' separately, though this will also be done in another publication. 



4. Numerical Results 

Before proceeding to show our results, we should mention the parameters we set for our 
simulations. We have taken the top (anti)quark to have a mass rrit = 171.6 GeV. The Z 
mass used was Mz = 91.18 GeV and was related to the mass, Mw, via the SM formula 
Mw = Mz cos Ow, where sin^ Ow = 0.222478. The Z width was F^ = 2.5 GeV. Also notice 
that, where relevant, Higgs contributions are included with Mh = 115 GeV. For the strong 
coupling constant, as, we have used the two-loop expression with Aq^^'*^ = 0.325 GeV in 
the MS scheme, yielding af^{M^) = 0.118. 

All the numerical results presented in this section are obtained considering a realistic 
experimental setup. Namely, partonic momenta are clustered into jets according to the 



Cambridge jet algorithm |36| (e.g., when yij < i/cut with i/cut = 0.001), the jets are required 
to lie in the central detector region 30° < 0jets < 150° and we request that the invariant 
mass of the jet system M^j is larger than 0.75 x ^/s. If a real photon is present in the 
final state, it is clustered according to the same algorithm, but we require that at least 
three "hadronic" jets are left at the end (i.e., events in which the photon is resolved are 
rejected). Notice that this procedure serves a twofold purpose. On the one hand, from the 
experimental viewpoint, a resolved (energetic and isolated) single photon is never treated as 
a jet. On the other hand, from a theoretical viewpoint, this enables us to remove divergent 
contributions appearing whenever an unresolved real gluon is produced via an IR emission, 
as we are not computing here ©(osOem) one-loop QCD corrections to e"'"e~ — > qqj. Finally, 
we sum over the final-state quarks. In order to show the behaviour of the corrections we 
are calculating, other than scanning in the collider energy, we have considered here the 
three discrete values of -y/i = Mz GeV, ^/s = 350 GeV and -y/s = 1 TeV. 

In Fig. 0, the relative effects^ on the cross section (integrated within the experimen- 
tal cuts defined above) induced by different contributions to the order asa"^^ correction 
are plotted as a function of the CM energy, in the range from 150 GeV to 1 TeV. The 
curves represent the effect of the QED (virtual and real) corrections only, the effect of 
the gauge bosons self-energy corrections (given by the graphs of Fig. |l|), the effect of the 
non-factorisable graphs of Figs. ^ (d)-(f) and ^ with WW exchange^*^, the effect of the 
weak corrections with the non-factorizing WW graphs removed (labelled as "full weak - 
non-fact WW graphs") and the total effect as the sum of the previous ones: the total 
effect is increasingly negative, reaching the —13% level at 1 TeV. It is worth mentioning 
that, as far as the non-factorisable WW corrections (represented by the graphs in Figs. ^ 
(d)-(f) and P) are concerned, in the case of d, s and b final-state quarks, only the direct 
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The effects are relative to the lowest-order cross section. 
"This is a gauge invariant subset of the complete correction. 
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Figure 7: Relative effect on ttie integrated cross section due to different contributions to the order 
a = Qem correction, as a function of the CM energy. For setup and input parameters, see the text. 



diagrams (Fig. ^ (d)-(e) and Fig. ^ (a)) are present due to charge conservation, while, for 
u and c quarks, only crossed diagrams are present, if the sum over initial- and final-state 
helicities is taken. In the case of ZZ exchange, all the graphs survive, giving rise to a 
cancellation at the leading-log level between direct and crossed diagrams, which does not 
occurr for WW exchange. Hence, the big negative correction is due to the presence of 
the WW non-factorisable graphs, which develop the aforementioned large Sudakov double 
logarithms in the high energy regime. 

We then show the impact of the EW corrections on some differential distributions 
of phenomenological interest. We organize the plots by showing the tree-level contribu- 
tions and the higher order corrections in three different contributions: the purely weak- 
interaction contribution (labelled "weak 0(a)"), purely weak plus QED corrections, which 
are dominated by the above-mentioned ISR (labelled "exact ©(a)"), and the weak plus 
electromagnetic correction in which the leading logarithms have been summed (labelled 
"exact 0{a) + h.o.LL"). The figures show in the upper panel the absolute distributions 
and in the lower panel the relative differences with respect to the tree-level rates. In 
Fig. ^, the azimuthal angle of the leading jet^^ distribution is shown at the Z peak. Here, 
the largest contribution to the total correction (—30% or so) comes from QED ISR because 
of the radiative return phenomenon. The purely weak corrections are negative and at the 
level of 3 — 4%. The higher-order QED radiation tends to compensate the order oem effect, 
since it enhances the cross section by ~ 5%. By raising the CM energy, the relative weight 
of the weak corrections becomes more important (see Figs. ^ and IC for the leading jet 
angle distribution at -y/s = 350 GeV and 1 TeV, respectively) and the effect of higher-order 



'^The leading jet is defined as the most energetic one. 
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Figure 8: Leading jet angle distribution at the Z peak. 
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Figure 9: Leading jet angle distribution at 350 GeV. 



QED corrections becomes negligible. It is worth noticing that, far from the Z peak, the cut 
M^j > 0.75 X y/s is more effective in reducing the radiative return phenomenon, reducing 
in turn the relative effect of ISR. 

In Figs. 11, 12 and 13 the leading jet energy is presented at ^/s = Mz, 350 GeV and 
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Figure 11: Leading jet energy distribution at the Z peak. 



1 TeV. This distribution is much more sensitive to the real radiation, as it can be expected. 
The corrections are very large in the distribution tail, where however the cross section is 
small. 

In the following some event shape variables are considered: the thrust T |37], the 
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Figure 12: Leading jet energy distribution at 350 GeV. 
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Figure 13: Leading jet energy distribution at 1 TeV. 



spherocity S [38|, the C-parameter [39| and the oblateness O [^]. One should notice that, 
owning to the fact that - as intimated - we are not computing here 0{asa'^^) one-loop 
QCD corrections to qq'y final states necessary to remove singularities associated with IR 
real gluon emission from qqg'j ones, we ought to maintain a i/cut between jets (we use 
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Figure 14: (1 — T)^ distribution at tlie Z peak. 
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Figure 15: (1 - r)ff distribution at 350 GeV. 

our default value), thereby defining the shape variables using the jets three-momenta, 
which at the same time enables us to remove energetic and isolated 'photon jets' from the 
final state. (In fact, for consistency, 'photon jets' are not used in the calculation of the 
shape variables even when the gluon is resolved in a qqg'j final state.) This is not unlike 
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Figure 16: (1 - T)^ distribution at 1 TeV. 



typical experimental procedures, whereby resolved photons from either the ISR or from 
the final state (i.e., those not originating from parton fragmentation) are separated from 
the hadronic system by using a jet clustering algorithm^^. In Figs. |l4[ ^ and ^ the 
distribution (1 — T)^ is shown. T is defined as 

T = max„^^J^, (4.1) 

\Pi\ 

where the sum runs over the reconstructed jet momenta and n is the thrust axis. The 
T distribution is one of the key observables used for the measurement of as in e'^e~ 



collisions [42|. It is worth noticing that while the purely weak corrections give an almost 
constant effect on the whole T range, the presence of the real bremsstrahlung gives a non 
trivial effect in the region T > 0.92. In view of a precise measurement of as at a future 
LC, EW corrections can play an important role. 

In Figs. 17, 18 and |l^, the distribution ^ is shown. S is defined as 



S = [-] mm^^^— — — , (4.2) 

V^y 22i\Pi\ 

where the sum runs over the reconstructed jet momenta and n is the spherocity axis. The 
effect of pure weak corrections is flat over the whole S range increasing in size from few 
percent at the Z peak to about 10% at 1 TeV. The addition of QED radiation at 0{a) 



^^This procedure is in fact standard at LEP2, see, e.g., Ref. ^ij, where resolved pliotons from both the 
initial and final states are more numerous than at LEPl, where the Z width naturally suppresses the former 
and where the phase space is more limited for the latter. 
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Figure 18: ^ distribution at 350 GeV. 



and at all orders does not change the picture except for the Z peak, where the 0{a) 
contribution is about 30% and the higher order corrections give an additional 10% effect. 
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Figure 19: ^ distribution at 1 TeV. 
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Figure 20: distribution at the Z peak. 
In Figs. 20, 21 and 22, the distribution ^ is shown. C is defined as 



^ ^ 3 Ei,j[\Pi\\Pj\ - jPi-Pj) /\Pi\\Pj\] 

"2 (EMf 



(4.3) 



where the sum runs over the reconstructed jet momenta. The effect of pure weak corrections 



- 21 - 



Born 
weak 0{a) 
exact 0{a) 
exact 0(a) + h.o. LL 




0.3 0.4 0.5 

C parameter 



Figure 21: ^ distribution at 350 GeV. 
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Figure 22: distribution at 1 TeV. 



is similar in magnitude and shape to the case of spherocity and thrust. The addition of 
QED radiation is very large only at the Z peak (~ 30%) while for higher energies it 
introduces a non trivial shape in the region 
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Figure 23: ^ distribution at the Z peak. 



o 

b 



12 
10 

8 
6 

4 ^ 

2 



Bom 
weak 0{a) 
exact 0{a) 
exact 0(a) + h.o. LL 







-1 




-| 




-5 - 






-10 - 


b 






b 


03 


-15 - 


1 

b 


b 


-20 - 






-25 - 




0.2 0.3 

O (oblateness) 



Figure 24: ^ distribution at 350 GeV. 
In Figs. 23, 24 and 25, the distribution ^ is shown. O is defined as 



O — F — F ■ 
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Figure 25: distribution at 1 TeV. 
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(4.4) 



where nmajor is an axis which lyes in the plane perpendicular to the thrust axis and the 
magnitude of the momenta projected into this plane is maximal along this axis; nminor is 
orthogonal to the thrust axis and nmajor- The numerical impact of the various contributions 
is similar to the case of the C-parameter. The region where QED real radiation introduces 
non trivial effects is C ^ 0.1. 



In Figs. 26, 27 and 28 the y distribution is shown. For each event, the observable 
y is defined as the minimum (Cambridge) yij such as yij > ycut = 0.001. Also on this 
distribution the weak effects are quite constant, while real radiation distorts its shape. 
Finally, Figs. 29, 3C and 31 represent the cross sections integrated over y in the range 
Vcut < y < Vmax, as & function of ymax- Corrections can be very large in both distributions 
generally at any energy, reaching the several tens of percent. 



5. Conclusions 

In summary, we have shown the phenomenological relevance that the calculation up to 
©(asa^Yv) '^^^ have in the study of (unflavoured) three-jet samples in e^e~ annihilation, 
for all energies ranging from ^/s = Mz to 1 TeV. Not only inclusive jet rates are affected, but 
also more exclusive distributions, both global (like the event shape variables) and individual 
(like energy and angles) ones. Effects range from a few percent to several tens of percent. 
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Figure 26: Cambridge y distribution at tlie Z peak. 
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Figure 27: Cambridge y distribution at 350 GeV. 



depending on the energy and the observable being studied, and we have shown cases where 
such higher-order contributions would impinge on the experimental measurements of jet 
quantities. 

A careful analysis of actual three-jet data is in order then, which should not only 
involve two-loop QCD contributions but also one-loop EW ones. One caveat should be 
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Figure 28: Cambridge y distribution at 1 TeV. 
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Figure 29: Integrated cross section as a function of the maximum y and relative corrections at 
the Z peak in the Cambridge scheme. 



borne in mind though concerning the latter. That is, as emphasised in Ref. |43|, particular 
care should be devoted to the treatment of real and Z production and decay in the 
definition of the jet data sample, as this will determine whether tree-level and Z 
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Figure 30: Integrated cross section as a function of the maximum y and relative corrections at 
350 GeV in the Cambridge scheme. 




Figure 31: Integrated cross section as a function of the maximum y and relative corrections at 1 
TeV in the Cambridge scheme. 



bremsstrahlung effects (neglected here) have to be included in the theoretical predictions 
through 0{asa'^yi), which might counterbalance in part the effects due to one-loop 



and Z virtual exchange. However, given the cleanhness of leptonic jet data samples, as 
compared to hadronic ones, we believe that the former contribution can always be effectively 
disentangled in data. 

We are now in the process of computing the ©(osckem) one-loop QCD corrections to 
e+e~ — > gg7, so that - upon combining these with the computations performed in this 
paper - we will eventually be in a position of rendering our complete predictions through 
that order of immediate experimental relevance. 
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Appendix 

In this Appendix we display the analytic expressions for the IR divergent part of the relevant 
scalar integrals used in this work. Such divergences occur either in the case in which all 
the particles inside the loop are taken to be massless (zero- internal-mass integrals), such 
as box or pentagon graphs in which both the exchanged gauge bosons are photons, or 
graphs which contain one massive internal particle (one-internal-mass integrals) but are 
nevertheless IR divergent, such as box or pentagon graphs with one photon and one Z 
boson exchange. We perform the integrals in (4 — 2e) dimensions with subtraction scale /x. 



The zero- internal-mass integrals have been discussed in detail in, e.g., Ref. |44|, but 
for completeness we reproduce the results here, introducing some compact notation for 
integrals that will be used to describe the one-internal-mass integrals: 



where 



Co(s, 0,0, 0,0,0) =— / 



p?=pi=o, {pi+P2)2 



r^(i-e)r(i + e) 

r(l - 2e) 



(V)"^ = r(l + e)(V)"Ml-^'^r + 



(A.l) 



(A.2) 



lf°)(s,t) = C7o(s,i,0,0,0,0) 



vr^ 7 k'^{k + piY{k + pi+p2)'^ 



pj=t, p2=o, {pi+P2)2 
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e^s-t) 



(A.3) 



= I)o(s,t,p', 0,0, 0,0, 0,0,0) = ^ 



r(e,/x) 

p2 



St 



2 \ ' 



-2Li.^l--J-2Li.^l--J-ln^-J --|. (A.4) 
For the above box integral, pf = p^, pf = 0{i > 1), (pi +^2)^ = •s, (p2 +^3)^ = t- 



4°'°^ (•512, S23, S34, S45, S51) = £^o('Si2, S23, S34, S45, S51, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) 



7r2 J k'^{k+piY{k+pi+p2Y{k-pi-p^Y{k-p^Y 



Pi=0, Sij=2pi-pj 



%^f^^^^VV2Li2(fl-^V2Li^(fl--l-? 

S23'S34S45 € V ■S5lSl2/i'' / V ■^Sl / V ^12/ 6 



+ cyclic permutations. 



(A.5) 



For the one- internal-mass- integrals, we begin with the IR divergent triangle integrals, 
i.e., those in which the internal massive propagator is adjacent to (as opposed to opposite) 
an external "mass". 



4^'^^(s,M2) _ C7q(s,o,0,M2,0,0) 

1 r 



(A;2 - M2)(A; + pi)'^{k + pi + P2f 



M2 

7^ 



1^'^ ((s-M2) 



(A.6) 



')(s,t,M2) = Co(s,t,0,0,0,0) 



r(e,/i) I J_ I f hP -s 



1 

^ 7 (A;2 - M2)(A;+pi)2(A;+pi +p2)2 
M2 - t 



Li2 



+ Li2 



(s - M2) 



(t - M2) 



(A.7) 



We now show how the IR divergent part of an n-point integral can be obtained in terms 
of (n — l)-point integrals. This method is similar to that of Ref. |45], but the resulting 
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expressions are somewhat simpler. Consider an n-point integral with general incoming 
external momenta pi (i = 1- ■ ■ n) and a massive propagator on the line jo, 



IniPi ■■■Pn,---'m 



■30 



where 



Qj = ^Pi, 10 = 0, 



(Ai 



(A.9) 



i=l 



and at least rrij^ is non-zero (there may, in general, be other internal finite masses, rrij). 

After Feynman parameterisation the integral over the loop-momentum k may be per- 
formed yielding 



^(Pi • • -Pn, • • -"ijo • • •) = r(n + e-2) 



n-2+e ■ 



(A.IO) 

The integration over the Feynman parameters can only give rise to an IR divergence when 
OLjQ = 0. In other words, the integral 



JO 



n-2+e 



(A.ll) 



is finite. 

The term Uj^ in the numerator can be generated by considering the integral 

1 f d(^-^'^k2G-^]qrk 



where the Gram matrix, G, is defined by 



Gij — 2qi ■ Qj. 



Now, by using the "pinch" relations 

2qi- k = [{k + qi)"^ - mf) - {k^ - ml) + mf - ml - qf 



(A.12) 



(A.13) 



(A.14) 



and the fact that the integral ( [A.12 ) is IR finite, we obtain the result for the IR divergence 
part of the integral I„, 



^hl ('^^ -ml- qf) I /„(pi • • - Pn, 



|IR 



-I (An _r{0} 
jol V (n-l) \n^l)J ' 



(A.15) 

where ij^^-i) ~ l)-point integral obtained by "pinching out" the Z-th propagator 

from the n-point integral. In- 

Some of the integrals on the RHS of Eq. ( A.15| ) may be IR convergent and may therefore 
be dropped. This process may then be iterated so that eventually the finite part of any 
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integral with one internal mass can be expressed in terms of integrals with zero internal 
masses and the IR divergent triangle integrals. 

The IR convergent contributions can then be obtained numerically, by taking the differ- 
ence between the original IR divergent integrals on the RHS of Eq. ( A.15| ) and regularising 



X 



the divergences by inserting a small mass for all the internal propagators and then check- 
ing that the result is insensitive to the value taken for this regulator mass, provided it is 
much smaller than any of the other masses or momenta but not so small that it introduces 
numerical instabilities. 

For the integrals required in this calculation, we have two IR box integrals with one 
non-zero internal mass and one non-zero external mass, depending on whether the external 
mass is opposite (a) or adjacent (b) to the massive propagator. 



1 r d^^-^')k 



7r2 J k'^{k+piY{{k+pi +P2Y - M'^){k - p4^f 
with pI=p^, pf = 0(i > 1), ipi+P2f = s, {p2 +P3)^ = t. 

lf\s,t,p',M^) = Do(s,t,0,/, 0,0,0, 0,m2,0) 

1 f d(4-2e)^ 



(A.16) 



(A.17) 



7r2 J k'^{k + piY{{k +pi+ P2y - M'^){k - p^Y 

with p\ = 0, = p^ pf = (i > 2), {pi+ p2f = s, {p2 + pzf = t. 

Using the above technique the IR divergent parts of these integrals are given by 

Aa). . 2 ^2. _ sli'''\s,M') + it-p')4'^'\t,p') 

I, {s,t,p ,M )iR - (A.18) 

and 

/<"(M.P=.M=).K = ^4"'(^.M-) + *4-'°W-pVf''(.p-.M-) , 

Finally we have the pentagon integral with no external masses and one internal mass, 
chosen arbitrarily to be the line between incoming momenta P2 and p^, 

li°'^\si2, S23, S3i, S51, M'^) = £^o(si2,S23,S34,S45,S51,0,0,0,0,0,0,0,M^,0,0) = 

p?=0, Sij=2pi-pj 



7r2 J k^{k + pif{{k + pi+ P2Y - M^)k -p^- pr>f{k - p^f 
The IR divergent part of this integral is found to be 
-f5'\si2, S23, S34, S45, S51, M^)iR 



S51'S45(S14S25 + 512^45 + 2si2Sl4 " S51S24 " 2M'^ S^) 
|(s5iS24 - S14S51S25 + Sl4S51'S24 " S14S25 " S12S51S45 + S12S14S45) 4"^ ('^IZ) ^23, S45, M^) 

+ (SI2S5IS45 - S51S24 + SI4S5IS25) lf\si2: S^l, Szi, M"^) + 2514551545/4^'°^ (551,545,523) 
+ (S51S24S45 + S14S25S45 - S14551524 + 514^25 " S12S45 " ■S12S14545) 4"^(s34, 523, ^51, M^) 
+ (si2545 — S51524S45 — 514S25545 — 2514551S45) /4'''*(534, 545, 5l2, M2)}. (A.21) 
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Erratum to: 



Full One-loop Electro- Weak Corrections to 

Three-jet Observables at the Z Pole and Beyond 

[JHEP 0903:047,2009] 

CM. Carloni-Calame^'^, S. Moretti^, F. Piccinini^ and D.A. Ross^ 

^ INFN, via E. Fermi 40, Frascati, Italy 

^ School of Physics and Astronomy, University of Southampton Highfield, Southampton 
SO 17 IB J, UK 

^ INFN - Sezione di Pavia, Via Bassi 6, 27100 Pavia, Italy 

All the integrated and differential cross-sections presented in the original paper have been 
mistakenly multiplied by an overall factor of 1/3. Furthermore, besides the cuts described in 
section an additional cut was imposed to obtain the results which is not correctly stated 
in the paper: namely, the energy of each jet is required to be larger than min(5 GeV, 0.05 x 

However, the relative effects due to the impact of the full one-loop EW terms of 
©(asOEivj) remain unchanged as well as the conclusions of our study. 
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